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ABSTRACT 


A techinque is described for the design of a multivariable 
controller for a VTOL (Vertical Takeoff and Landing) helicopter, based 
on a recent technique incorporating pole assignment and minimum eigen- 
value sensitivity. A fourth order state space model of the helicopter 
available in the literature and used by earlier researchers is used in 
the present design also. Since the elements of the system matrix of 
the model do change over varying flight conditions, the minimum sen- 
Sitivity design approach used in this thesis has a decided advantage 
over many previous design procedures which neglected the parameter var- 
jations and designed fixed feedback controllers or resorted to more 


complex adaptive techniques. 
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perturbation in 


approaches or takes the value of 


characteristic equation coefficient 


less than 

greater than 
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frequency 
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used for matrix identification 


used for matrix identification 
derivative 
pole(eigenvalue) location 
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CHAPTER (1) 
INTRODUCTION 


1.1 History of Vertical Take-off and Landing (VTOL) Aircraft 

One of the first concepts of VTOL* aircraft was recorded by 
Leonardo da Vinci in 1483. It was not until near the end of the 
eighteenth century when the first physical model of a VTOL aircraft was 
made by Launoy and Bienvena [15] .»* 

During the next one hundred and twenty years there were many 
concepts proposed and implemented on numerous VTOL aircraft configurations 
in many countries of the world. In 1907 Breguet of France designed and 
built a VTOL rotary wing aircraft which was successful in lifting itself 
and a pilot off the ground. Although this aircraft could fly, it could 
do so only when tethered to the ground since the pilot had absolutely no 
directional control over the aircraft. It was this lack of control which 
plagued every inventor and engineer who pioneered in the development of 
the VTOL aircraft. 

1.2 Development of Control Systems 

During the first World War, Vonkarman and Petroczy produced a 
capative helicopter used only for army observation duties. This helicopter 
could be stabilized only by use of its mooring cables. From 1910 to 1920 


the advances seen in the development of the VTOL aircraft was almost nil. 


* Hereafter in this thesis, VTOL will be understood to mean specifically 
helicopter. 


*k Numbers in rectangular brackets refer to referances listed under 
Bibliography. 
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The next ten to fifteen years saw many experimental designs of 
helicopters. During this period of time the designers used a "direct" 
method for pilot control over the various modes of operation of the air- 
craft. This is shovn in figure 1 for only one control mode, forward 
flight. Although this was a feasible method of control, it was not prac- 
tical over even a short period of time because constant and absolute 
pilot attention was required. 

In 1935-1936 Breguent-Dorand of France was able to maintain a 
helicopter airborne for areasonable time which exhibited good control 
characteristics. Th: Breguet-Dorand helicopter provided, 

i) cyclic pitch control - to govern horizontal flight; 
ii) collective pitch control - to govern vertical flight; and 

iii) differential collective control - to govern yaw motion. 
Rather than the direct control, these controls were connected by systems 
of bell-cranks, levers and plates from the pilot controls to the respec- 
tive controlled mode of operation. This type of control was used and 
improved upon over the subsequent years. 

It was not until the introduction of modern control theory that 
a significant change was seen in the design of controls for the VTOL air- 
Cratius 
1.3 Application of Modern Control Theory to VTOL Aircraft Controls 

With the introduction of modern control theory, the control of 
the VTOL aircraft was greatly increased and improved upon. This was pos- 
sible since now the pilot was able to divert his attention to tasks other 
than strictly that of flying the aircraft. 


In the modern day helicopter the input commands of the pilot 
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Figure 1. Direct Control 


are supplemented by means of feedback from the output of the control 
system. This feedback can be accomplished electronically in a number 
of ways. 

A constant gain feedback controller, which is nearly always 
preferred because of its simplicity and cost, is one means of providing 
the necessary control and stability augmentation. With the increasing 
advancements in technology more sophisticated control systems are being 
proposed. For instance adaptive controllers are used, which in many cases, 
require on-board computers for continual updating of the feedback matrix 
to ensure stability of the system [1] . These systems while being very 
promising, are also very expensive to implement on existing VTOL aircraft. 
Until the recently proposed systems can be implemented at a more resonable 
cost there may still be ways of re-designing the existing systems so that 
they are more acceptable to the aircraft industry both from the point of 
view of performance and cost. 

What follows in subsequent chapters outlines a new procedure 
for designing a simpler and less costly constant gain feedback controller. 
The method results in a more stable system with better response charac- 
teristics over the flight regime than those of some previous design pro- 
cedures which also employ constant gain feedback controllers. 

1.4 Scope of Thesis 

The aim of research reported in this thesis is to show the 
feasibility of using pole assignment and minimum sensitivity design 
technique to develop a controller for use in a VTOL aircraft. 

In Chapter 2, a brief outline of the method to be used will be 


given. Also an arbitrary second-order linear time invarient multivariable 
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system will be used to demonstrate the approach used. The results will 
then be compared to those obtained using an aribitrary feedback controller. 

Chapter 3 shows the results obtained using the minimum 
sensitivity approach when applied to a helicopter. In this study the 
model used by K.S. Narendra and S.S. Tripathi [2] is used. An arbitrary 
controller is then used on the same system and the results compared to 
those obtained using the minimum sensitivity design approach. 

Chapter 4 discusses the introduction of integral control 
along with minimum sensitivity design in developing a controller for 
improving the transient and steady state performance of the VTOL air- 
Cratt. 

In Chapter 5, the results are summarized and conclusions 
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CHAPTER (2) 


POLE ASSIGNMENT WITH MINIMUN EIGENVALUE SENSITIVITY 
TO PLANT PARAMETER VARIATIONS 


2.1 Introduction 

In this chapter a procedure is described for assigning the 
closed loop poles of a feedback control system to specified locations in 
the complex frequency plane in such a way that the poles have minimum 
sensitivity to plant parameter variations. The discussion follows closely 
the work of Gourishanker and Ramar [24] . This design procedure is used 
later in this thesis for designing feedback controllers for VTOL aircraft. 

In the design of control systems, mathematical models of the 
plant or process to be controlled are used. In many cases the parameters 
of the plant are subject to variations caused by several factors such as 
environmental changes, external control inputs, ageing of components, etc. 
In many cases, even small changes in the values of the plant parameters may 
affect the system behaviour appreciably. Since the design of controllers 
js usually based on nominal values of plant parameters it becomes necessary 
to study the effect of parameter variations on system behaviour. The 
sensitivity of system behaviour to variations in plant parameters as re- 
flected in the mathematical model has engaged the attention of reasearchers 
in recent years. One measure of this is known as eigenvalue sensitivity 
which is defined as the sensitivity of the closed-loop poles of a system 
to variations in plant parameters. Although eigenvalue sensitivity is a 
less direct measure of system performance it is being recognized more anc 
more as a very useful measure of system performance [18, 24, 26] . One 


advantage of using eigenvalue sensitivity is that it can be easily 
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computed. It is well known that all the closed-loop poles can be 
assigned exactly to desired locations when all the state variables are 
available for feedback [2g] . However, the solution is not unique in that 
more than one feedback controller can be designed to achieve the same 
objective as far as pole locations are concerned. Different feedback 
controllers may result in different behaviour of the system from some 
other point of view. This design freedom can be used to satisfy some 
other design criterion without altering the desired pole locations. 

In the procedure described in this chapter this design freedom 
is used to minimize the sensitivity of the assigned closed-loop poles 
with respect to variations in the parameters of the plant. In other words, 
pole assignment and insensitive design are combined. The closed-loop 
poles are specified a priori in order to achieve a certain transient re- 
sponse and the feedback controller is designed to achieve these pole lo- 
cations. The sensitivity of these poles to parameter variations will be 
a minimum. 

2.2. Design of a Unity Rank Feedback Controller for Pole Assignment 

Consider a linear time invarient multivariable system repre- 


sented by 


where x is an n-dimensional state vector, u is an r-dimensional input 
vector, A is the n x n system matrix and B is the nxr input matrix. This 
multi-input system is reduced to an equivalent Single input system by 


defining a new scalar input u' as 


where 
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q =f], lye eee yg GQ] wee nen n nee eee + ee (2-3) 


is an r-dimensional column vector, and R is the r-dimensional external 


input control vector. Then using (2-2) equation (2-1) becomes 

3 . 

Keer NX ee BOU ee B Ree eee ene ee cee eee eee (2-4) 
Now a feedback control 
is designed such that the poles of the closed-loop system 


x = (A + Bqk)x + BR 


fe) 


x= (A+ ba k)x + BR )---------------------------- (2-6) 
spe aes 
x = A + BR 
where 
b = B 
a . 
A=A+b-k 
q 
are at the specified locations S12 Soo +++ Sp. Note that k is an 


n-dimensional row vector 

k= [kyy Kyo cee Kyy] -eenerteeeccccceseseeeeee (2-7) 
and is a function of the r-dimensional vector q. In general q is chosen 
arbitrarily. In this thesis q wil be chosen to minimize the sensitivity 
of the eigenvalues of the closed-loop system. Once q = chosen the vector 
k is determined for the specified pole locations. The rxn unity-rank 
feedback matrix K for the given closed-loop system (2-6) is easily 


obtained as 
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2.3 Eigenvalue Sensitivity 


Since we are using unity-rank feedback, the vector q will be 
chosen to ensure minimum sensitivity of the eigenvalues of the closed- 
Toop system, thus resulting in the minimum sensitivity feedback controller 
as defined by (2-8). 

Following Morgan [26] the sensitivity oH of any pole 


(eigenvalue) s. of the closed-loop system matrix A in equation (2-6) with 


‘ 
respect to a "small" variation in the element aso of the open-loop system 


matrix A is defined as 


where 
g'(s.) is the derivative of the closed-loop system 
characteristic polynominal with respect to 
Ss evaluated at s = S33 


RUS. ese the adjoint(s.1 - A), I is the nxn 
Identity matrix; and 


tr denotes the trace of a matrix. 


As can be seen from (2-6), the closed-loop system matrix A, is 


a function of q and k;_ k is a function of q and therefore the sensitivity 


j 


Sao is directly a function of q. In order to minimize the effect of the 
variations in the parameter aso on a a sensitivity functional is form- 
ulated [26]as 
eo i ee 
ey DAUR Qi eee pee gre (2-10) 
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J obviously is a function of q. The objective is to determine 
q such that J is a minimum. Note that (2-10) is general in that it is 
apdlicable when all the elements of A are subject to variation. Now to 
facilitate the minimization of J with respect to q, the sensitivity 
functional, J should be evaluted in terms of q. This is accomplished if 
the sensitivity, Si can be determined as only a function of q. One 
method of achieving this is to transform the matrix pair (A,b,) to 
phase-variable or canonical form. By using a method reported in the 
literature on a transformation matrix P(q) is determined which trans- 
forms the matrix pair (A,b,) to phase-variable form, i.e the trans form- 


ation 


OK i a a a mae ee me ee oe ee meee ee ees em oes me ome a (2-11) 


is obtained. 
Now using the transformation equation (2-11) one obtains for 


the transformed closed-looped system equation 


2 = P(q)(A + b.k)P(q)7! z + P(q)BR 
or to eee (2-12) 
Z = Aoz + P(q)BR 
where 
ho = P(q)(A + byk)P(q) © =2---na===-—-a---n- (2-13) 
The preselection of the desired poles of the closed-loop system matrix A 


determines the elements of Aco, where Ao in canonical form is written as 
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where 1s %» ++. a are the known coefficients of the closed-loop 
characteristic equation 


1 n-2 


Gi) tua tS S bl gp eee On S + Oy) = QO --- (2-15) 


n n-1 


We can rewrite (2-9) in terms of the transformed equation as 


where 
se pS) .tne adjoint(s.1 - Ao), noting this is 
independent of q by the definition of Ao 
in (2-14). 


Now substituting (2-13) into (2-16) we have 
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Since only the elements of A are varying 
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then using equation (2-18) in (2-17), one obtains for the sensitivity 


where 


is the %th row of P(q) 7! 
) is the jth column of P(q) 


functional, J strictly as a function of the variable q. We can now 
minimize the functional J with respect to q; the required k necessary 
for the desired pole locations can then be easily determined by use 
of equation (2-13). Once k is known, the minimum sensitivity controller 
K is obtained by equation (2-8). 
The above procedure is utilized throughout the remainder of 
this thesis to obtain the minimum sensitivity controller. 
2.3 An Example 
In this section a numerical example is presented to illustrate 
the minimum sensitivity design procedure described in the previous section. 
Consider a linear time-invarient second-order system with two 
inputs and two outputs described by 


(@) 


X = AX + Bu ---------------- - +--+ eee (2-20) 
where 
1 a 1 
A = ; BS 
i 1 1 
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Now using the unity-rank feedback matrix as defined by (2-5) and (2-8), 


equation (2-20) becomes 


x = (A+ bgk)x +8R 
aeheere CUT! 9 § ~ eee eee (2221) 


Ax +BR 


x 
Since there are only two inputs we can write q in the following 
form 
q= [1 gy) ---nnnnnnnnnnnnonnnoee- (2-22) 
where qy the scalar is to be determined such that the sensitivity 
functional J is a minimum. 
It can be shown that the matrix which transforms the pair 
(A,b,) is the same as that which transforms the matrix pair (A ,b,) to 
phase-variable or canonical form [36] . 
To determine the transformation matrix, P(q), the controll- 
ability matrix Q is formed where 


q= [apt ab) =e (2-23) 


or 


We then have 
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The trarsformation matrix 


is obtained as follows: 


i) set the last row of ie =p 


1 
ii) solve the equation Po = PyA. 
Therefore 
atest W 
WS | eereners - 
q il 2q4 + 1 
Po = a eee 
2 aay ae 2q, +1 
Se pa ce 
eq, ap wh 2q4 ual 
P(q) = 
“apts! ven Lee 
2@q, + 1 2q, #1 
ard 
i q4 
Pia) a 
=| q + ] 


The poles of the open-loop system are 9 and 2 as deiermined 
by the roots of det(sI - A) = 0, which results in the open-loop equation 
S(s - 2) = 0 --------------------------------- (2-24) 
Since the open-loop system is unstable, we wish to stabilize 


the system by using a constant gain feedback controller which will be 
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designed using pole assignment and minimum eigenvalue sensitivity. In 
order to obtain a "good" transient: response we will choose new closed- 
loop system pole locations as Sy * -1 and Sy = -2. The characteristic 
equation of the closed-loop system is then 


(Sate) (Se eee =a 


? See eR nega eer ig errr ae ee ee ee tore (2-25) 
Sua sO Sauer) 
Comparing coef“icients of equation (2-25) and (2-15), we get 
a, = 2 
TG ip) eee eet eer, Zhe ona Bea (2-26) 
Oy = 3 
In canonical form the closed-loop system matrix then has the form 
0 il 0 1 
Ao = = | | wae nee eeen--- (2-27) 
“Oy “Oy -2 -3 
Now 
S; ia 2) 1 
Se) Baad US leew ic a ety 
i j 
-2 Sa 
7 
fOlaaige= 4 eg -1, we get 
2 1 
L(s,) = | | ww ee ene eee nnn -------------- (2-28) 
-2 -l 
and for i = 2; oe Sy = -2, we get 


Now g'(s.) is the derivative of the closed-loop characteristic equation 


(2-25) evaluated at Se (ij = 1,2), therefore 
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and for i = 2 


(2-19) 
Sea (akeuier ac) ee ee (2-33) 
ues g'(s.) Q ieee 
TOY wen 12 
a he 
2 = 1,2 


then for i = 1, using the results of P(q), P(q)7!, (2-28) and (2-31) 


_ (2qr + 3)(q, - 1) (qq + 3) (ay + 2) 


2q4 at 294 Tel 
at 2 
A 
(20a eel (2qi0- tq) = 2) 
2q, 2q, +1 


and for i = 2, using the results of P(q), P(q)7!, (2-29) and (2-32) 


Next we shall consider two cases of plant parameter 


variations. These are manifested as variations of the elements of the 


system matrix A. 
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Case l 


Nenote the matrix A as 


aa 499 
and assume the element Qo4 alone is subjected to variations. Then the 


sensitivity functional to be minimized is 


ea 
tt 
iW f&-INoO 
— 
— 
=A 
— 
PO 


(Say 
H 
— 
v2) 
— 
~— 
DO 
+ 
— 
NY 
NO 


By inspection the minimum value of J = 0 occurs if q, = SO la eintinite 
values of q, ore disreguarded for practical reasons. Choosing q, = 1 


gives the follwing results: 


: 1 
q= [1 1] aly 
ae 
Si get at 
P = : cae 
eer 


Wr 


i 
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Now it is known that by using the transformation matrix P, the canonical 
form of A is Ao. Therefore using the relationship shown in equation (2-13) 


we can solve for the only unknown value, k, which assigns the poles to 
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the desired locations in the conplex frequency plane. Thus we have 
Nc ON ere cerca cece (2-35) 
and using equation (2-8) we find the minimum sensitivity controller kK 


to be 


This feedback matrix, K, should result in the poles of the closed-loop 
system to be at the desired locations of s, = -l, S5 = =e IS can 
easily be varified by taking det(sI - (A + BK)) = 9. 


Now let the elenent G5, vary such that hay, = 0.2, then 


) 0 
Na 
OFZ ) 
Denote 
Aaa eo aoe erent een ae ec ee (2-37) 
or 
i 1 i 
A = 
Ce ] 
Then the perturbed closed- oop system matrix is 
Tcl Nae CMe Mme ee A act (2-38) 
or 
~ -? 0 
A = 
-48 -1 


which results in a closed-!oop characteristic equation of 
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of ergs “hae y= 0) 
~------------------------- (2-39) 
(Sere?) (Site 1) =U 
indicating we have closed-loop pole locations of Sy = -1 and So = -2. 


Thus we see there are no variations in the pole locations when the 
element Ay varies from 1 to 1.2. In fact, as can be seen from fi any 
variation in A545 will have "no" affect on the pole locations. 

Now suppose we disregard the sensitivity criterion and assign 
the closed-loop poles to the same locations as before by choosing an 
arbitrary q. Let q = [1 ones Then following the procedure 


described earlier we get 


27S: 

ae a) 
- 270 150 

fal 21 


By taking det [s1 - (A + BK) = 0, it is easily shown that the 
poles of the closed-loop system using arbitrary feedback are at the exact 
locations. 

Now using the same variation in An, as before one obtains for 
the closed-loop characteristic equation 

pera Pr (07 14 oe I yemeceen ea te, Jaen ene eeeeee (2-40) 
which has pole locations of S) = -0.13 and S5 = -2.87. Thus we see that 


for a 20% change in An4> Sy changes by 87% and So changes by 43.5%.* 


* See Appendix III for the method used to calculate the percent change 
in pole locations. 
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Case 2 
Now suppose all the elements of A are subjected to variations. 
The functional to be minimized is then 


n n n 


SPO ae dslape 
i=l j=l e=1 J* 
1S USA Reo CSalr 2 ea 
Za we CC TR YS 2ei\e 
Sich Una Aba leh ce elas 8: 
and using the values obtained for Si and s the above becomes 
2 4 3 2 
hte r oi ee + 50) -- (2-41 
(2a, * 1) % yet teas 32q, + 76q5 + 48q, + 50) (2-41) 


Equation (2-41) was solved on the digital computer at the 
University of Alberta center using the minimization Subroutine ZXMIN 
(Int. Mathematical Statistical Lib., Vol. 1, Ed. 5) to determine the 
minimum valve of J. The minimum value of J = 47.097 occurs when 
q, = 0.714. Then using the value of q = [1 0.714] | the minimum 


sensitivity controller was determined as 
-3.471 -1.472 
si i -1.05 

This feedback matrix obtained using the minimum sensitivity design assigns 


the poles of the closed-loop system to the exact locations of S, = -1 and 


Taking the variation in A to be 
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results in th: perturbed closed-loop systen matrix of 
ny -1.378 -.15 
A = 

-4,849 -1.422 


which has pole locations of 


Thus we see that Sy changes by 45% and So changes by 12.5%. 
Now using an arbitrary controller by choosing gq = {i 10]' 


one obtains (as before) the feedback matrix of 


= al 15 
Clee et 
Kos 
Pec) ue toe 
al 2h 


Then subjecting A to the same viriations as shown in (2-42) and using 


the above K, one obtains the po’es of the perturbed closed-loop system as 


0.14 
= -2.94 


ai 
a) 

indicating that Sy changes by 114% and S. changes by 47%. Also as can 
be seen by the above pole locations, the arbitrary system becomes unstable 
when A is subjected to the variztion indicated. 

Thus one can see that the minimum sensitivity controller design 
has definite advantaces over the arbitrary feedback controller design 
for pole assicnment. 

Both systems were simulated with the external control input 
Peon io ‘bon for the cases AS = 0 and AA = equation(2-42). The 
responses obtained wren using the minimum sensitivity controller are 


Shown ir figure 2 for xy and in figure 4 for Xo3 when uSing the arbitrary 


controller the responses are shown in figure 3 for xy and figure 5 for Xo. 
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Figure 2. 


WITH AA=O 
—e—e~ WITH AAO 


eee = 
TIME (sec) 


Response of X, Using the 
Minimum cence? Controller 


a 


—— WITH AA=O 
—e-e WITH AAO 


ey aibis ois 
TIME (sec) 


Figure 3. Response of X1 Using the 
Arbitrary Controller 


ra) 


IRE fA SO) 
a WITH AA#O 


8 
TIME (sec) 


Figure 4. Response of Xo Using the 
Minimum Sensitivity Controller 
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| 
—— WITH AA=O 
-—e WITH AAO 
Py) 


| vemetaiober 


R 4 6 8 lO 


| TIME (sec) 
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Figure 5. Response of x9 Using the 
Arbitrary Controller 


CHAPTER (3) 


VTOL CONTROLLER [ESIGN 


In general, the dynamics of a VTOL aircraft are such that it 
is inherently unstable. For this aircraft to be of any practical use, 
some form of stability and/or control augmentation is required so that 
it will be stable and controllable during its flight regime. The intro- 
duction of feedback controllers into the VTOL system is one way of 
achieving this. 

Over the past fifty years or so many controllers have been 
designed and proven to be useful on existing VTOL aircraft. Some of 
these are outlined in the references at the end of this thesis. 

Although stability is perhaps the main criterion in the design 
of a VTIOL aircraft contro!ler it is not by far the only factor to be con- 
sidered. Response of the system to predetermined or disturbance inputs, 
cost of implementing the control system and reliability of the system 
are also important criteria to be considered when designing a controller. 

It is shown in this thesis (in this chapter and the next) that 
a constant gain feedback controller designed by using the minimum eigen-~ 
value or pole assignment enhances the stability and transient response 
characteristics of the VTOL aircraft considered, even under conditions of 
plant parameter variations. Also the cost of implementing a constant 
feedback control scheme (as presented here) is relatively small compared 
to some of the other more sophisticated control schemes, such as the 


adaptive controi scheme proposed by Narendra and Tripathi [1] ‘ 
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3.2 Model Dynamics 


The model used in the design procedure is that used by 
K.S. Narendra and S.S. Tripathi [1] . As in most other techniques used 
for VTOL controller design this model is obtained by linearization of the 
system dynamics around a nominal air speed. 
The linearized model of the VTOL aircraft in the vertical plane 
is described by 
; a) UWE eS) SEE Se Ge ee eee ec (3-1) 
where A is the (4x4) system matrix, B is the (4x2) control matrix, x is 
the 4-dimensional state vector and u is the 2-dimensional control vector. 
The state variables are: 
Xi 7 horizontal velocity; 
Xo - vertical velocity; 


pitch rate, and 


x< 
(Se) 
1 


Xq > pitch angle. 
The control inputs are: 
Uy - collective, and 


Uy - longitudinal cyclic. 


* The control u, is located on the collective pitch lever at the pilots 
SiGe hes matn use is the control over the vertical velocity of the 
VTOL aircraft by the selection of a desired flight angle. The flight 
path or angle is changed merely by the pilot applying collective 
control pressure in an up or down motion. This control also has some 
affect on the horizontal velocity. 


The control u, is one of the controls located on the cyclic control 
stick immediately in front of the pilot. Its main use is to control 
the horizontal velocity of the VTOL aircraft. The up and down movement 
of the control (the one we are considering) changes the forward hori- 
zontal velocity of the helicopter while a left or right movement of 
the control will change the heading reference to that of the flight 


indicator. 
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In this thesis, following [1] , the nominal airspeed is assumed 
to be one hundred and thirty-five (135) knots. At this nominal airspeed 


the system matrix A and control matrix B are: 


122 943 344 -0.0366 0.0271 0.0188 -0.4555 
0.0482 -1.01 0.0024 -4.0208 
D7 LOCZ 0s 3601s Oly Oy, 1.42 
0.0 0.0 is) 0.0 


0.4422 0.1761 
3.5446 -7.5922 
-5.52 4.49 
0.0 0.0 


As the aircraft speed deviates from the nominal air speed all 
the elements in the first three rows of both matrices change. The most 
Significant changes occuring in the elements Ao and azn» the rest of the 
elements can be assumed to remain constant without serious loss of 
accuracy .* 

The design technique described in chapter (2) is used here to 
take into account the changes in the parameters Azo and Azq as the air- 


speed changes from sixty (60) to one hundred and seventy (170) knots. 


* In this thesis the design technique is used only to take into account 
the changes in the elements of the plant matrix A. However, if found 
necessary, this design can easily be modified to incorporate the 
changes in the control matrix B as well fig} . 
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It is assumed that the elements A35 and Az, vary in a linear manner as 
the airspeed changes from sixty to «ne hundred and thirty-five knots, 
and from one hundred and thirty-five to one hundred and seventy knots. 
Although the above assumption is not entirely valid, it is acceptable 
to illustrate the design procedure. The manner in which 830 and asa 
vary as a function of time is shown in figures 6 and 7 respectively. 
Note that the time shown in figures 6 and 7 is only approximate as it 
had to be taken from the graphs of reference 1, Page 194. 

Since two elements of the matrix A are changing with changes 
in airspeed the functional to be minimized is 

d= COE Wp RNY “acest a (3-2) 
42 

Since the closed-loop system is of 4th order, there are four eigenvalues. 


The desired locations for these eigenvalues are 


Sy = -1.5 
So = -2,0 
S = AP Sean 
Sq = -] - jl 


The Fortran program listed in Appendix I was used to carry out 
necessary computations for the entire design procedure. This program was 
written to incorporate the Subroutine ZXMIN which computes the minimum of 
the junctional J (3-2). 

The minimum value of | = 0.0054 was obtained when 1. O05 6/7; 


This resulted in the minimum sansitivity feedback controller 
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Figure 6. 


TIME (sec) 


Change in the Parameter a, of the 
System Matrix A as a function of Time 
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-3.268 -9.1334 0.8567 2.930 
K = --- (3-3) 


-2.212 -0.0903 0.5799 1.984 
At the nominal airspeed of one hundred and thirty-five knots, this feed- 
back controller (3-3), assigns the closed-loop poles exactly to the desired 
locations. 
Now using the minimum sensitivity controller (3-3) the system 
was cheked for stability and pole variation of the closed-loop system 


at the extremes of the flight reqime speeds, i.e. at sixty and at 


one hundred and seventy knots. /it sixty knots Ago = 0.06635, a3, = Oak he) 
and at one hundred and seventy knots Azo = 0.5047 and zy = 20208 
At sixty Knots: 
Denote 
0 0 0 0 
0 0 0 0 
V0) | no Nn) beer t ooo (3-4) 
0 Mazo 8) Aaa 
0 0 0 0 
where 
Aa», = -(),30175 
hag, = -1.3002 
Let 
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AV) 
Now replacing A by A in equation (2-6) one obtains for the closed-loop 


system 


" 
(A + BK)x + BR 


Pa) XS) Ss 1S) 
il 


Ax + BR 
2 

where A is the closed-loop system matrix, and R is the external control 

vector, the magnitude of which is determined by the pilot. 


~% 
When the pole locations of A were checked it was found that 


$1 = -0.796 
Sy = -1.55 
S2 = SALAS SE sean basis 
S, = =158 sj. 23 


which indicates that the system is stable at the lower extreme of the 


flight regime. Thus one can see that for a Aa,, = -0.30175 which 


32 
corresponds to a 82% change in a5 and for a hag, = -1.3002 which 


corresponds to a 92% change in An the closed-loop poles change by 


AS, = 0.704 or a 47% change in Sy 

AS5 = 0.45 or a 22.5% change in S5 

AS, = -0.58 + j0.33 or a 46% change in s, @ /4.9° 
AS, = a25S- 10.33 0% de4o, changes in Sa Q {-4.9° 


At one hundred and seventy knots: 
Using the same notation as before 


Aago = ORG 


hazy = 1.106 


alr 
Again checking the pole locations of A the closed-loop system resulted in 


pole locations of 
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Sy = -1,52 
Sy = -2,56 
33 = -0.71 + j0.9:16 
Sq = -0.71 - j0.9¢6 


Thus for a change in a0 of 37% and a change in An of 78&% the 


closed-loop pcles changed by 


AS, = -0.02 or a 1.3% change in Sy 
AS, = -0,56 or a 28% change in So 
AS, = 0.29 - j0.054 or a 20.8% change ins, © /-8° 
bS_ = 0.29 + 50.054 or a 20.8% change ins, © /& 


Again we see that the system is stable at tie upper end of the VTOL flight 
regime. 

By uSing the standard method for pole assignment incorporating 
unity rank feedback one can obtain for the VTOL system by choosing arbi- 
trarily, a= [1 10] !, the feedback controller 

12063) ©-O 9302.8 ==(es09 -0.535 
10,63. =3. 02 ~7 09 ~5.35 

At the nominal air speed, it can easily be shown the arbitrary 
controller of (3-7) assigns the closec-loop poles to the exact desired 
locations as in the previous instance. 

AL Sixty KnOLs:. 


The c!osed-loop poles were found to be 


Sy = -2.92 + j4.65 
Sy = -2.92 - j4.65 
S3 = Q. 16/7 4° 30.764 


0.167 - j0.764 
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From the location of the poles of the closed-loop system it is clear 
that the system is unstable when the plant matrix A is subjected to the 


variation of Aa3z5 = ~0.30175 and haz, = -~]1.3002 which resulted in pole 


location changes of 


= 
2) 
li 


j = 71.42 + 54.65 or a 121% change in s, @ /-57.9° 


AS, = -0,.92 - j4.65 or a 118.5% change in s, ®@ /579° 


> 
iva) 
! 


3 = 1.167 - j0.236 or a 84% change in s, @ /-57.3° 


AS, = 1.167 + j0.236 or a 84% change in s, @ / aay 


for a change in a5 and a4 of 82% and 92% respectively. 
At one hundred and seventy knots: 


Tne closed-loop poles were found to be 


Sy = -1.026 
Sota -5.83 
S2 = O41 (oa sg 24, 
Jae Omi iG wices 


Again by the location of the poles in the complex frequency plane one can 
see that the system is unstable when the plant matrix A is subjected to 
the variation of hazy = 0.1365 and hax, = ',106 which resulted in pole 


location changes of 


AS, = 1.474 or a 98.3% change in Sy 
AS, = -3.83 or a 91.5% change in So 
as, = 1.176 + j1.4 or a 129% change in s, @ £49.2° 
AS, = 1.176 - j1.4 or a 129% change in s, @ /49,2° 


for a change in Ao and a, of 37% an'l 78% respectively. 


Thus one can see that the design procedure which minimizes 


eigenvalue sensitivity has definite advantages over an arbitrary controller 
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which also assigns the closed-loop poles to the desired locations at the 
nominal speed of operation. 
3,0) (racking Capabilities of the Controller 

The ability of a control system to alter the steady state 
condition of the aircraft to a new steady state condition in a satis- 
factory manner, as determined by pilot control inputs, is a measure of 
how well the system is able to track or follow desired command inputs. 
An over damped system response will result in the aircraft being too 
Sluggish or slow to reach its new desired steady state condition where as 
an underdamped system may result in uncontrollable oscillations of the 
aircraft when it is subjected to command inputs. A schematic diagram of 
the VTOL control system in the tracking mode is shown in figure 8. 

By choosing the dominant poles of the closed-loop system as 

S30 7 -i+ jl 

at the nominal airspeed of one hundred and thirty-five knots results in a 
damping ratio of ¢ = €.707 and the undamped natural frequency of oscilla- 
tion a = 1.414 radians/second. At sixty knots the damping ratio x = 0.765 


and o, = 2.065 radians/second, and at one hundred and seventy knots ¢ = 0.6 


if 


and Wp Lelee radians/second. 


The above values of cz and w, are within the boundries of 
a Ole Ne ore h 2 critical and 0.5 radians/second < pies 5 radians/second as 


determined by Buffum and Robertson [10] . 


To check the tracking capabilities of the minimum sensitivity 


controller and the arbitrary controller the following procedure was 


adopted, 
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As mentioned previously, the model being used is that used by 
Narendra and Tripathi [i]. From [1] it is impossible to determine the 
exact pilot input which corresponds to any particular air speed. For 
this reason, in the following simulations, the pilot inputs (in this case 
R) were choosen as step inputs representing a particular horizontal 
velocity (see Table I). The pilot inputs are used to Jjenerate the 
corresponding control input to the system which results in an arbitrary 
steady state condition for each of the four output states X19 Xo X3 and 
Xn (the arbitrary steady states are obtained since we are arbitrarily 
selecting the pilot input to represent a varticular air speed). 

At sixty knots when Azo = 0.96635 and a 


ge 5 
input was chosen to be R = [0.4 0.4] i at one hundred and thirty-five 


es ).1196, the pilot 


knots Azo = 0.3681 and Az, = 1.42 the pilot input was taken as 
R = fees 0) y, and at one hundred and seventy knots Azo = 0.5047, 
Az, = 2.526 and the pilot input was taken as R = [1.25 1.25} °, 

The system was simulated on the digital computer using the 
CSMP prcgram listed in Appendix II to determine the steady state values of 
the variables corresponding to each of the three steady state conditions 
mentioned in the previous paragraph. 

It was now assumed that the aircra*t was flying in steady 

state at sixty knots with R = [0.4 0.4]! and initial conditions on the 
variables as determined previously. Also by prior simulation at one 
hundred and thirty-five knots the steady state values of the variables at 
this air speed were also determined. 


The pilot input: was now suddenly increased from R = [0.4 0.4]! 
Cones [1.0 1.0]! to simulate an acceleration from sixty to one hundred 
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and thirty-five knots. Since the parameters Azo and 234 change with 
changes in air speed, the change in the system matrix A was programed 

into the simulation (this was done by modifying the program in Appendix II). 
The response of the aircraft was then determined and plotted as shown in 
figures 9, 10, 11 and 12 for X42 Xoo Xq and Xa respectively. 

The above procedure was then repeated to simulate an acceler- 
ation from one hundred and thirty-five to one hundred and seventy knots. 
The response was again obtained and plotted as shown in figures 13, 14, 

15 and 16 for Xp Xos X3 and Xq respectively. 

Now for comparision purposes an arbitrary feed matrix was used 
(the one determined in Section 3.2). The previously determined steady 
state values of X19 Xoo Xq and Xa at sixty, one hundred and thirty-five 
and one hundred and seventy knots respectively were used to determine the 
respective pilot input required to achieve the same steady state values at 
the desired air speeds when the arbitrary controller is used on the same 
aircraft. It is felt this is valid reasoning since for a completely 
different controller the pilot inputs will be different to achieve the 
Same steady state results. Using this method it was found that to obtain 
the same steady state values as previous at sixty, one hundred and thirty- 
five and one hundred and seventy knots that R = [-0.59 Se 
R= [21.25 -13.9]' and R= [21.45  -15.15)' respectively. 

The system was then simulated starting in steady state at 
sixty knots (R = [-0.59 -7.68]") and suddenly accelerating to one 
hundred and thirty-five knots (R changed to R = [-1.25 Higeg |e The 
response was obtained, plotted and compared to the response obtained under 


the same conditions using the minimum sensitivity controller shown in 
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Figure 9. Response of Xj when Accelerating from 
Sixty to One Hundred and Thirty-five Knots 
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Figure 10. Response of X9 when Accelerating from 
: Sixty to One Hundred and Thirty-five Knots 
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Figure 11. Response of X3 when Accelerating from 
: Sixty to One Hundred and Thirty-five Knots 
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The above procedure was then repeated to simulate an accel- 
eration from one hundred and thirty-five to one hundred and seventy knots 
using the appropriate pilot input and initial conditions. The results 
obtained were then plotted and compared to those obtained for the same 
Situation when using the minimum sensitivity controller, see figure 13, 
14, 15 and 16. 

It must be noted that although we are concerned in this design 
mainly with controlling the horizontal speed of the aircraft that the 
control inputs used, namely collective and longitudinal cyclic, produce 
changes in the other variables as well as the horizontal velocity 
variable X1- This is common in all VTOL type aircraft as the modes of 
operation of the VTIOL aircraft are not sufficiently decoupled so that 
changes in one mode of operation will not affect changes in another [2]. 
The longitudinal cyclic input from the control column is the main control 
over the horizontal speed of the aircraft. The collective input from the 
collective pitch lever affects not only the horizontal velocity but also 
the vertical velocity. Just how one affects the other is quite compli- 
cated and beyond the scope of this thesis. The references listed at the 
end of this thesis provide ample information on this aspect [15] . 

3.4 Regulating Capabilities of the Controller 

The regulating ability of any controller is quite important, 

especially in a VTOL aircraft. By regulation is meant the ability of the 


controller to maintain a steady state condition in the presence of dis- 


turbances. 
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Figure 13. Response of Xj when Accelerating from 
: One Hundred and Thirty-five to One Hundred 
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Figure 14. Response of Xg when Accelerating from 
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Figure 16. Response of Xq when Accelerating from 
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From equation (3-1) which describes the dynamics of the heli- 
copter in the vertical plane we have designed a unity-rank minimum sen- 
sitivity controller of the form shown in equation (2-2). For regulation 
purposes we will redefine (2-2) to take into account disturbances by 


defing 


where 


u' = equation (2-5) 


Rae plot. input 
Ww = appropriately dimensioned disturbance vector 

In all previous calculations and simulations w was taken as 
identically zero (whether w = 0 or w # 0 has no effect on the design of 
the controller). It is now desired to see how the controller will re- 
gulate the speed of the VTOL aircraft in the presence of disturbances, 
inevew 720), 

Assume that the pilot is flying the helicopter along in a 
steady state condition at one hundred and thirty-five knots. Suddenly 
there is a gust of wind such that w = Cw, Wo = [0.4 0.4)! which 
lasts for three seconds. Without any corrective pilot action the response 
of the system variables of the closed-loop system with the minimum sensi- 
tivity controller shown in figures 17, 18, 19 and 20 for X12 Xos X3 and Xa 
respectively. 

For comparison purposes the closed-loop control system with an 
arbitrary controller having the same pole locations as the minimum sen- 


sitivity controller was subjected to the same test. The results are shown 


a 


Vit 7 


ae i; a 


Seah ae 
ee 


on - La 
is methine © SUERTE i 


ere saa) naan mo 


iamig lille 79 Inieee RP? ae wh 


Pe oN ' ‘hap , e 


‘ con Be Haws 
7 ae 7 
= ° i - _ i. i 7 
(Put) antcnupe 0°@ AF) at aoe 
: 2) 
Kat & i > 
jae? dee © 7 » i 
- 
oo9 iets Te anol VT od op TQeree ¥ 7 P i 
3 
Lwig ba catty ahh ween iter. Ghee th a 
Ye mplashoed 4 Co +i sin eagle 


cas She ~etfl aa i bey pe jc cavteee or a? Amtieam a Bt 
On. cae med. .4) IDK ry a any " waa 


a r a q 
, “ae tp ei 1 - q v; 


fp onris sodas v@ OAT UNIS! af al gf tile 

ingen a? sda <¢t t+ ys ial perrintse shi Ay, sale AE 
aL ees ie dl = ere ou shoe ee 
ipmaaens pls. tation Jeti ay ttatiren via ‘poeta POI 
irene mantis git tt Ye wees sea -obedlng 9 , te : 
pein oh sak te) HS Wilt At jeer al # 


i 


MINIMUM 
SENSITIVITY CONTROLLER 


-~@—e ARBITRARY CONTROLLER 


“15 
TIME (sec) 


lO 


Disturbance Effect on X1 when Flying in 
Steady State of One Hundred and Thirty-five 


Knots 


Figure 1/7. 


51 


nenigetiaga NS 8.7 
AB JSGR7 NOD VR ART ah 7 


& wien CO ie ne irs 


MINIMUM 
SENSITIVITY CONTROLLER 


—@-—-e ARBITRARY CONTROLLER 


TIME (sec) 


i 18. Disturbance Effect on X9 when Flying in 
a Steady State of One One Hundred and Thirty-five 


Knots 


52 


7 
¥ 
= 7 
| 
= 


| = Eee : 
5 
_ 6 
% 7 - _ 78 7 - 


J - 
: 7 = 
ia 


; sh San : 
a) ie ey, . 
ae a et y 


2.0 


2.0 


-1.0} 


“1.5 | 


MINIMUM SENSITIVITY 
CONTROLLER 


~®—e- ARBITRARY CONTROLLER 


10) Se eners 
TIME (sec) 


oa 19. Disturbance Effect on X3 when Flying in — 
gue Steady State of One Hundred and Thirty-five 


Knots 


53 


MINIMUM | 
SENSITIVITY CONTROLLER 


~@—@ ARBITRARY CONTROLLER 


TIME (sec) 


Figure 20. Disturbance Effect on Xq when Flying in ; 
Steady State of One Hundred and Thirty-five 
Knots 


54 


MEIIGATHOS: 7H eG = ae: i. ; 
ne usbeYHOS esi SS | x 
a 


ee ee ee 


Be il aD. ON BS 
o) ‘“! 


Tesh aah 


aye 


in figures 17, 18, 19 and 20 for X12 Xo» X3 and x, respectively, to make 
comparison easier. 

From the results shown in figures 17, 18, 19 and 20 it can be 
seen that the regulating capabilities of the minimum sensitivity controller 


are far superior to those of the arbitrary controller. 
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CHAPTER (4) 


INCORPORATION OF INTEGRAL CONTROL 


4.1 Introduction 

In the previous chapter the controller was designed without 
any consideration given to the steady-state errors produced when the 
system is subjected to step command inputs. Steady-state errors will 
always occur when a proportional control system is activated by step 
command inputs. One way of reducing the steady-state errors is to in- 
crease the loop gains considerably, which in general, is undersirable and 
unacceptable as far as the present work is concerned since the controller 
is being designed for minimum eigenvalue sensitivity. The steady-state 
error can also be eliminated by the inclusion of integral control action 
in the controller. This method is preferable here since integral control 
can be incorporated in the design procedure. 

When the problem was first tackled, it was hoped that the 
steady-state error of all the four state variables of the helicopter 
system under study could be made zero. However the present state of the 
art does not allow the augmentation of the original system by more than 
the number of inputs to the original system [ie] . The system under con- 
sideration is of fourth order with two command inputs. Hence the integral 
control can be applied to not more than two of the four states of the 
system. In view of this restriction it was decided to apply the integral 


control to the horizontal (x,) and the vertical (x,) velocity of the 


aircraft system. 
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4.2 Minimum Sensitivity Controller Design With Zero Steady-state Error 


Consider a linear multivariable multi-input dynamic system 


described by 


where 

xX = n-dimensional state vector; 

u= m-dimensional input vector; 

y = p-dimensional output vector; and 

A, B.and C are respectively nxn, nxm and pxn matrices. 

It is now desired to augment the original system such that the 

transient response of the controlled variable(s) is/are satisfactory 
(response characteristics to satisfy some previous design specifications) 
and that in the steady-state the output variables become equal to some 
arbitrary constant reference inputs, i.e , there is no steady-stare error. 


Using integral control we introduce new state variables 


where 
z= p-dimensional vector; 
X= p-dimensional state vector to be controlled;and 
I= pxp identity matrix. 


Combining equations (4-1) and (4-2) the augmented system is 
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For the single input case it is a fairly simple matter to 
choose a feedback matrix that will assign the poles of the closed-loop 
system to the desired locations if the system given by equation (4-4) is 
controllable. 

For a system of the form (4-4) with more than one input, the 
system matrix A* will never be cyclic. since it has multiple eigenvalues 
of zero. Thus it is not possible to choose a unit-rank feedback matrix 
such that the poles of the closed-loop system will be at the desired lo- 
cations in the complex frequency plane. Thus the design must proceed in 
two stages. 

First a feedback matrix omar is chosen so that the first 
stage closed-loop system matrix is cyclic. In order to facilitate this 


we choose a control law of the form 
ue = Te + u2 + Ro -nnnn nnn nnn nnn ene ener eee (4-5) 


where ul and ie are respectively the first and second stage controls and 


R is the pilot control input(R is an mxl input vector). Define the first 


** The matrix A* will be cyclic if there exists an n-vector w such that 
the matrix [w A*w ... A*"AJis of full rank. 

**k Note that Kt, can be randonly chosen but KI must be of full rank p. 
The choice of Kt, can be used to satisfy other design specifications 
not pertaining to minimum sensitivity analysis. 
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stage control as 


where 
ere AOU id eer eae es (4-7) 


and when p = m, as in our case, a first stage control matrix of the form 
shown in (4-7) will generally result in the removal of the non-cyclicity 


of the system matrix if 


king Kiyo oe Kip 1 Om: 0 
KI = Kin, King ait Kion = 0 1 . O| ----- (4-8) 
ki kin... ki 
40 uy, Ki np 0 0 1 


and KI is pxp. 
Now using equations (4-4), (4-5) and (4-6) results in the first 


stage closed-loop system equation 


4 
x A '-BKI} |x] [Bl |B 
at eee ore a ol seh tes | arta] Reet ace x 
S Ei On| |b 0 0 
ee (4-9) 
- 2 
x = A*x + BuO + BFR + eh 


where the definition if A* is obvious and B* and I* were previously defined. 
The control scheme is shown in schematic form in figure 21. 

Now the first stage closed-loop system matrix, Ax, of (4-9) can 
be considered as the open-loop system matrix for the second stage design. 
The design procedure used in chapter 3 can now be used to assign the poles 
of the closed-loop augmented system using unity-rank feedback and minimum 


sensitivity analysis by choosing a second stage control of the form 
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where q is the same as in chapter 2 and 


1S Vee Beatie ral)! 9 sess RS ee rte etn (4-11) 


and Kyo is an n+p row vector defined as 


Ke oee. | kere emkd ete kn me are KPinep | -777- (4-12) 
which we shall denote as 
{ 
Keo = Lkd | kp] 0 --------------- 2-22 enone nee ~-- (4-13) 


Haine utd Te -ckester 10 tne spocme SiS do. Je meneepnon (4-14) 
' tk 
ul i= -Kro [x : rl 


where Kr5 is the minimum sensitivity controller defined by 


Bagh hs Fas... a el eee eee omc (4-15) 


K £2 


TZ 
The total control Jaw is then obtained by use of equations 


(4-5), (4-6) and (4-14) and defined as 


or 


where 


Recalling the definition of be in chapter 2 and noting the definition of 


above, we obtain the total closed-loop augmented system 


Kee anaan 


we T2 


equation 
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fe) ' ; 
X A - Bae PBK I= bekp) 1x B 0 
; ae Zhen 0 : (4-17) 
mene 
xX = AX + BYR + leae 
where the closed-loop system matrix is 
' 
. |A - b-kd 1-BKI = b kp 
Wie | aye tea) aa en Le || ee (4-18) 


-C 0 
The closed-loop poles of the closed-loop system matrix A can 
be placed in the complex frequency plane to any desired locations by the 
procedure used in chapter 3 if the system of (4-9) jis controllable. 
The system is controllable if and only if: 
i) the controllability matrix 


Q= Lb, : AMD! ae yr | Be eae (4-19) 
is of full rank n; and 


ii) the matrix AX is of rank n+ p. 
For proof of this see Young and Willems [16] . 
4.3 Integral Control For VTOL Model 
The procedure described in the previous section will now be 
applied to the VTOL system under study. The horizontal and vertical 
velocity of the VTOL aircraft are chosen as the output variables to 
which the integral control is to be applied. 


The two new variables used to augment the original system are 


then described as 
(@) 
Se wa | 


ak 


ul 


Ven r2 
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where 
hy We horizontal velocity; 
Xy = vertical velocity; 
pe reference horizontal velocity; and 
eos reference vertical velocity. 
Stage 1. 


We choose a first stage control law as described by equation 


(4-6) by selecting a first stage control matrix Ky, as 


As in Chapter (3) we again design the controller at the nominal air speed 
of one hundred and thirty-five knots. With Kr4 as above the first stage 


closed-loop system matrix as defined in equation (4-10) is 


-0.0366 0.027 0.0188 -0.4555 -0.4422 -0.1761 


0.0482 -1.01 0.0024 -4.0208 -3.5446 7.5922 
“ 0.1002 Os36819 F-07707 1.42 boo -4.49 
stig 0.0 0.0 120 0.0 0.0 0.0 

-1.0 0.0 0.0 0.0 0.0 0.0 

0.0 -1.0 0.0 0.0 0.0 0.0 


We now consider Ax as the open-loop system matrix of the 


augmented helicopter system. 
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Second Stage. 


It is now desired to obtain the second stage controller using 
minimum sensitivity design to assign the closed-loop poles of the augmented 


helicopter system to the desired locations in the complex frequency plane 


of 


$1 = ~4.9 
So = -3.0 
$3 = -2.0 
San -1.5 
S¢ = -1.0 + jl.0 
open pects Ja. U 


Since only the parameters a5 and Az, are varying with changes 
in the horizontal air speed the minimum eigenvalue sensitivity functional 
is given by 


n * : 
ee CS ces.) ee (4-20) 


where the sensitivity Sly and oi were defined in Chapter (3). 

This functional is a function of q. The Fortram program 
listed in Appendix I was modified to accommodate the augmented sixth-order 
system, then used to do the entire design calculations. The minimum value 
of J occurred for q = [1 0.562]'. This then resulted in the minimum 


eigenvalue sensitivity feedback matrix Kyo of 


8.8193 OnOU28 Be ora loyg -10.1264 -9,6836 8.3666 


K = 
T2 49553 © 0.3381 1.2012 -5.6897.  -5.4408 4.7009 


* 


Pies on be 


: es : é oy 


@ - 
a ee - 


e 
ait sa ee : Ma 
ghey WAT ee roo aeaye saa = 

wet digce at YS whley tut bones ony, aah 


inte yin ae! wiges ait. at: oh at 


ein “a 


on! ' : ear 7 — 


im eapddni arin 
coe] pecan | 


rt 


a 


At eae 

el ae nies y | 

woe Re ee cae any es eoore _ 
bom? ws otavennte wliapy. ain hah ment 2 


ee ? ere pat sues 

awed Ph EAM: Lb 4a ee Ca ebayaall'd ri 
a ee ee adeun ener cege : gas tVGkats hint 
wird by wigigln by ailT emit aber pdtv pil we) eo 
sivstnist 0 o' bed turer wa ADR: pea ae +3 


or the total ninimum eigenvalue sensitivity feedback matrix K. of 


t 
Sygoon aeIG: 0.6028 =Ce13/o -10.1264 -8.6836 8.3666 
4.9553 0.3381 -1.2012 -5.6897 -5.4408 5, 7009 


using the Ky above resulted in the closed-loop system matrix as defined 


in equation (4-18) of 


-4.8091 -0.2986 Te 57, 5.0244 4.7980 -4.7036 
6.4084 “0.57/60 -1.5394 -11.3237 -10.5281 13.6260 
26.9309 2.1718 -7.1149 -28.9313 -23.5041 20.5868 


A = 
0.0 0.0 120 Oe) O50 0.0 
-1.0 0.0 0.0 0.0 0.0 0.0 
0.0 -1.0 0.0 0.0 0.0 0.0 


which has the closed-loop pole locations at the desired positions in the 
complex frequency plane. 

To check the sensitivity of the pole locations to variations of 
the elements a5 and a4 of the matrix A, the pole locations of the closed- 
loop system matrix A were computed at the extremes of the flight regime, 


i.e. at sixty and one hundred and seventy knots, where a35 = 0.06635, 


Az, = 0.1198 and Azo = 0250475 Aon = 2.526 respectively. 
At sixty knots: Aaj, = -0.30175 or an 82% change in a5 4 and 
baz, = -1.3002 or a 92% change in a3: The closed-loop poles were found 


to be 
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-0.7652 
=O. /G0GRtedleoc 
-0.7606 - j1.589 


which indicate pole location changes of 


il 


=O. 07455 J  Slsaoy 
= be 0/4 Lleol 3emor 
-0.065 or 
0.7348 or 
0.3494 + j0.589 or 
0.3494 - jo.589 or 
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a 38% change ins, @ /-20.4° 
a 62% change in s, @ /20.4° 
a 3.2% change in S4 

a 49% change in Sy 

a 48% change in sp @ /-19.4° 
a 48% change in s, @ /19.4° 


At one hundred and seventy knots: Mazo = 0.1366 indicating a 


change in a5 of 37% and dar, = 1.196 indicating a 78% change in ao4- 


The closec-lcop poles were found to be 


-5.03 

arAR Ma! ce AME shes) 
=2,159 - 31.6083 
-2.075 

=, 5307 +7.)0.6379 


-0.5367 - j0.6379 


which indicat’ pole location changes of 
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-1.03 or a 25.7% change in Sy 
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AS. = -0.159 - j1.683 or an 84% change in S, @/-7° 


ce 
nn 
I 


ne -0.575 Or a 38% change in Sa 


Ass = 0.4633 - 50.3721 or a 42% change in sp @ /-4.9° 


[Be 
7p) 
i 


g = 0.4633 + 50.3721 or a 42% change in s_ @ /4.9° 


For comparison an arbitrary controller was designed to assign 
the poles of the closed-loop system to the same locations as before. 
Ghcrcetofq= (1 33) | resulted in the arbitrary feedback controller of 

Se SHAS) -1.2647 “3.0333. =/ 1020 -15.3930 3.4701 
2 
ie 39 .9526 -3.7940 -10.0000 -21.3059 -46.1789 10.4103 


when Key was the same as in the minimum sensitivity design. The complete 


feedback controller is 


e375 -1.2647 -3.3333 -7.1020 - 14.3930 3.4701 
K = 
| 39.9526 -3.7940 -10.9000 -21.3059 -46.1789 11.4103 


The arbitrary controller above resulted in the closed-loop system matrix 


of 
-12.9613 1.2544 3.2538 6.4370 14.4967 -3.5438 
256.1712  -25.3318  -64.1044 -140.6057 -299.5822 74.3289 
-105.7743 10.4231 25.7931 57.8806 127.8941  -32.0772 
he 0.0 0.0 1.0 0.0 0.0 0.0 
=1s0 0.0 0.0 0.0 0.0 0.0 
0.0 She 0.0 0.0 0.0 0.0 


which has the desired pole locations. 


Again checking the pole locations of the system when using the 


arbitrary controller and subjecting A to the same variations as in the 


67 


4 
or 
Y 
i 


- he.% 


. 
be * 


“* tv ae 7 


oo 
m0 


9.6. : 


> Sunnis a ‘ *t “rest 


pelted Sur ne 91, et ani vo oti 


Te fui ee ee ee oad toate ee 


aor 
o_rge . GD Mabe wep ee Mtey. ‘i 
a 4 han ; ic Tan et 


ee ce 


onels yt i vag pee: ut 


ey, 
. daw “OP? cosids patee “ at aria 


ae fe (S05 Sit- PRE 
h cthwts/ Page a ni eR?’ @ be eae wh’ fou oe 


ef aol Lovina? toa ey’ 
; 7 
nyt. CUA ae ae 7 


’ a ; 8 7 7" 
oie ere Vent”. » Ue : 
AT eae jv omg? oy itt, nt ay Be: iv: " at Tear na od ar 


¥ =o | 
\¢@ Sade iv | 7 


ibe are tt Ao état. s-.] ; 
vot. ry 


~~ 
7 


minimum sensitivity case, it was found that at: 


Sixty knots 


Sia -9 326 
So = -1,85 
S3 = -0.4069 + j0.5479 


Sq = ~0.4069 - j0.5479 


Se = aJNpraejoey Go alsincists 
Se = -0.2555 - j3.386 
resulting in a 
AS, = -5.326 Or a 133% change in Sy 
AS, = Seek) or a 72% change in S5 


AS, = 1.5931 + j0.5479 or an 84% change in s, @ /-8.4° 
As, = 1.0931 - j0.5479 or an 82% change in Sa @ /8.4° 
Asp = 0.7445 + j2.386 or a 177% change in sp @ /-40.8° 


eee 


ASe = 0.7445 - j2.386 or a 177% change in sp @ /40.8° 


One hundred and seventy knots 


S, = 74.955 + 34.145 
Soy ah LEIS j4.145 
S$, = -1.936 

Sp = -0.4509 

ot aa (CoO e weet lel a7 
Sie -0.086 - j1.177 


resulting in a 


iis -0.955 + j4.145 or a 106% change in Sy @ /-40° 
ASy = 1.955 - j4.145 or a 153% change in s» @ /40° 
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AS. = 0.064 or a 3.2% change in S4 


AS, so deatouke. or a 70% change in Sa 


0) 
- @ /-40.8 


ASg = 0.914 - j0.177 or a 66% change in s_ @ /40.8° 


AS = 0.914 + j0.177 or a 66% change in s 


A comparison of the above results reveal that the variations in 
the poles of the closed-loop system, due to changes in the system plant 
parameters, are appreciably less when the minimum sensitivity design approach 


is used to obtain the feedback controller. 


4.4 Tracking Capabilities of the Zero Steady-state Error Controller 

To illustrate the tracking capabilities of the two controllers 
designed in section 4.3 the system was simulated on the digital computer 
using the procram listed in appendix II. The simulation here is slightly 
different than that of chapter 3. 

Now a prior knowledge of the command inputs must be known 
(16 37] i.e. X, must be known. Then in the steady-state the controlled 
outputs (the outputs to which the integral control action is applied, in 
this case Xy and Xo) should approach the values of Xaq and X09 respectively, 


i.e., the steady-state error should be zero. 


The responses of both systems were obtained when accelerating 
from sixty to one hundred and thirty-five knots and again when accelerating 
from one hundred and thirty-five to one hundred and seventy knots. The 
following steps were involved in the respective simulations: 

i) to begin the simulation at sixty knots 
a) R= (0.4 0.4] T vas arbitrarily chosen; 


ne a Aol XJ and X..5 equal to 1.0 and 0.0 respectively 
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the steady-state values are obtained for the 


other vaviablescas, xo =.17 0s xe c= Oe Otex earn 


1 2 3 


and Xq = -0.018. 
ii) to accelerate from sixty to one hundred and thirty-five knots 
ay R= (1.0 1.0] I was chosen; 
b) initial conditions were set as obtained in (i); and 


Co Xn] and X no equaled 3.0 and 0.0 respectively. 


iii) to accelerate from one hundred and thirty-five to one 
hundred and seventy knots 


a) R= (1.25 1.25] ' was chosen: 
b) initial conditions were set as obtained in (ii); and 


c) X nq and Xn9 equaled 5.0 and 0.0 respectively. 


A schematic diagram of the control system used for the above 
simulations is shown in figure 22. 

The responses of the horizontal (x,) and vertical (x5) velocities 
are shown in figure 23 for xy and figure 25 for Xo when accelerating from 
sixty to one hundred and thirty-five knots and in figure 24 for Xy and 
figure 26 for Xo when accelerating from one hundred and thirty-five to 
one hundred and seventy knots for both the minimum sensitivity and the 
arbitrary controllers. 

4.5 Regulating Capabilities of the Zero Steady-state Error Controllers 

In this section we look at the ability of the controllers to 
maintain a constant output, i.e. zero steady-state error of the chosen 


outputs, when the system is subjected to disturbances. 


The original equation (4-1) is rewritten here to include the 


disturbance vector w as 
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XGe] AeutBUGeiBye ee ee (4-21) 
where w is previously defined. 
Let x, u, w be the nominal values of the respective variables 
and A be the nominal value of the system matrix A. In terms of the 
deviations of the variables from their nominal values let 


fe) O 0 
MM ae oe GP Oe 


A> A+ 6A 
Then equation (4-21) in terms of (4-22) becomes 
i | 
x + 6X = (A + 6A)(x + 6x) + B(u + Su) + B(w + dw) ---- (4-23) 


or the perturbed equation in terms of the deviations of the variables is 


sx'= (A+ 6A)ox + Béu + Bow + (GA)x -eb-2---2----- (4-24) 
where 
ob oe hy % 
Om) *0, » 0m r6 — 1% 
§ A = X = ooo 
0 baz, 534 X3 
Cio mow FD & 
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sy 
x ~OX 
Rly <1 Ngan cero! «Fe eee eee ae eee ee (4-25) 
6X ~8X, 
(since Xn] and Xo are the reference values and therefore 6X14 = 8X9 = 0.) 


in order to maintain the deviation 6X4 and 5X5 of the horizontal and 


vertical velocity at zero in steady state. The augmented system then 


becomes 
Ms AUae % av A, ee 
6x = Aéx + Bou + Bow + (SA)X 9 ~---------~------ (4-26) 
where 
ey. 10 
a %v See ee ee 
6X = : hits c | 0 
B wean v |X 
B =|-- CLES ie in ces reed Mies. 
9) 0 @) 0 


Equation (4-23) represents the system simulated on the 
computer at the nominal airspeed of one hundred and thirty-five knots 
with. the perturbed control law éu of 


Uk = -K 8X own nnn enn nnn enn nen en ne nnnee (4-27) 


where K, is defined in section 4.2 and has the form 
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Since u = u_ + gu we have for the control law 


Guy = -kdy 46x, -kdyoSxq kd, 36x -kdy g6xXq -(kiy, + kP1q)oxet (kKipo + kPy 9) dX¢ 
SUy = ~kdy 16x, -kdyg6Xq -kdy,6Xq -kdygdXq -(king + kPoq)OxXG+ (king + kPo9) Xe 


The control scheme is shown in figure 27 where u is replaced 


by the arbitrarily chosen pilot input R where R is taken as 


For the simulation the nominal values of A and B were taken 
from Narendra and Tripathi [1] » and are listed in section 3.3. Fora 
desired horizontal velocity of one hundred and thirty-five knots and zero 
vertical velocity the following nominal values were selected 


i) pilot inputs 
Vis ihe 
ey iba 
ii) reference values 
Kay = 3.0 
Xo = 0.0 
For the desired horizontal velocity of x, = 3.0 and vertical 


velocity of Xo = 0.0 the following steady-state values were obtained for 


the other variables when: 
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i) using the minimum sensitivity controller 


a) X= 5 
Xn = Oe) 
X= O50 
Xq = -0.067 
b) u, (collective) = eile 


u, (longitudinal cyclic) = 0.116 


ii) using the arbitrary controller 


a) x1 = sEe 
Xo = On0 
Xz = 0.0 
Xq = -0.067 
b) u, (collective) = 0.132 


u, (longi tudinal cyclic) = 0.116 

Both systems were then simulated and subjected to the same 
disturbance vector éw = Low, bw] | as shown in figure 28. The responses 
of xy and Xv for both systems are shown in figure 29 and figure 30 
respectively. 

Upon comparison of the responses, X4 and Xs it is seen that 
the variations produced when using the minimum sensitivity controller are 
somewhat less than those produced when the arbitrary controller is used. 

It is also noted that when checking the responsed of xy and 
Xo in all the figures of this chapter, that because of the incorporation 
of integral control for the variables xy and Xo there is no steady-state 


error in any case (i.e., in steady-state the outputs x4 and Xo equal the 


reference inputs Xd and Xn9 respectively). 
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CHAPTER (5) 


SUMMARY AND CONCLUSIONS 


5.1 Discussion of Results 

A study of the results reveals that the minimum sensitivity 
design approach in constructing a constant gain feedback controller for a 
system has definite advantages over an arbitrary feedback controller when 
assigning the closed-loop poles of a system to specified locations in the 
complex frequency plane. 

Since the location of the closed-loop poles of any system is an 
important factor determining the stability of the system, the sensitivity 
of these pole locations to variations in elements of the plant matrix are 
also critical. The sensitivity design approach ensures minimum variation 
of the pole locations when the element(s) of the plant matrix vary from 
their nominal values. This is illustrates quite effectively in chapters 
Zee ones 

The poles of a closed-loop system are also a factor in deter- 
mining the response of the variables of a system. As can be seen by 
checking the responses in the previous three chapters, those obtained 
using the minimum sensitivity controller are much better than those using 
an arbitrary controller since the poles using the minimum sensitivity 
controller are more insensitive to change when the elements of the plant 
matrix vary, ie. they tend to remain at the desired locations. 

The results of chapter 4 further illustrate the effective use 
of the minimum sensitivity design approach. When using minimum sensitivity 


analysis, the «ffective incorporation of integral control to ensure zero 
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Steady-state error of the selected variables and also decoupling of the 
modes in steady state, would suggest that the minimum sensitivity analysis 
procedure is not restricted in its applications but that it can quite 
possibly be used in conjunction with a number of existing control methods. 
5.2 Possible Areas for Further Research 

The research in this thesis is intended to show one possible 

method of designing a constant gain feedback controller for VTOL aircraft. 
The results indicate that the controller can be effective over the flight 
regime of a VTOL aircraft. The following are some areas for further 
research. 

(1) Chapter 2, 3 and 4 took into account only the variations 
of the plant matrix A of the system. A further study 
could take into account the variations of the input 
matrix B and also variations in the feedback matrix K. 

(2) In this thesis it was assumed all the state variables 
were available for feedback, and in many cases this may 
not be the case. A study could be conducted for the 
case where some of the states are not available for 
feedback. The incorporation of an observer (state 
reconstructor) to estimate the unmeasurable state vari- 
ables could be considered. The minimum sensitivity design 
approach could then be used to design the feedback con- 
troller as before. 

(3) The use of some other method other than unity-rank feed- 


back could also be used to obtain the minimum sensitivity 


controller. 
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Fortran Program Used to do Design Calculations 
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APPENDIX II 


Continuous System Modeling Program (CSMP) Used For Simulations 


LABEL HELTEOPTER SYSTEM 

TNT T TAL. 

UR TeNOMTNAL PILOT COLLECTIVE INPUT, 
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X4O02 TNT AL CONCUITION FOR X4, 
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APPENDIX III 


Calculation of Percent(%) Change in Pole Locations 
Suppose the pole locations before variations in the system 


matrix are as shown in the complex frequency plane of figure 1. 
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Figure A-1. Pole Locations 


Before variation we will denote the poles as 


$1 = -5 
So = -4 
S = -2 
it's wl eraeae 
$s, = =] - jl 


Due to some change in the system matrix the pole locations 


change and the new pole positions are represented as 


191 


ru 


mefeys sad ai par ere wed protzenet ofa, oe neat? 
,Penetl? Fe actla @avutet aelenes: aft ot ars 20 oh! 1g don 


uf _~ ee 
a Le ‘ ° 
' er 
u ay : ! a 
~— ¢ Te 7 
i 
. * a) 7 
4 
> 
v7 a ee pr rx) 
» : 5 i ; i- be = 
' ? | ( , 
| mt | | 
1 
i i 
' Ham | ; 
i 5 
reed) alO8 f-N  67,.oi7 _ 
9 


> f 7 a ro - 
v titi, aA? wjaheat rit ak) nottnaay: te ia 


Sy = -5.5 

S5 = -3.4 + j2 
Spee co 4 geade 
sq> 1.0°+ 50-5 
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which indicate a change in the pole locations of 
AS, = -0.5 
AS, = OF6RE Ic 
As, = -1.4 - j2 
AS, = 200 
As, = Cmte Old 
The above changes are calculated by the equation 


ASen eS e=nS eel eae eee Whe eat acer 


j 1 j 


The percent change in Ss; is then calculated by 


% change of S; = alee coats 
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of poles 


Now if there is a change in the angle of the pole location 


as well, it was calculated as 
* 
ly JR DS i 5 


A positive anale was taken as shown below in figure 2. 
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Figure A-2. Positive Angle 


va 


<b == 9 n= Ls % *, sate 


Ayia + Ost ley “ye 


Bie Ra” 


ae - 


{ 2 4.4! @ 124 
3 . | P+ & i © a : 
| EERE mas 


; ‘ : a 
é fy 2 rh M4 © gfe < 


aed PabeseroN on eopnets.2 pie 


oh « wiiellw: > ei ee ST He y, - e feta 


0 ae) At aoe a # a. by ail eri ot 1 
ey 53 
~aties wehliety stale i tagee't v 7 


7 AL ay 


eK, * a4 + a ‘ 


A at a ei6y’ vw rqaeta 0 ation | . ite 


ria 
iy ri ° 
4 


» 
> 


eh rioh! ote a Te Siefia sieves pecs ay 


_ 7 


mt 
a « 


iy 


AS¢ 


O25 


x 100 


10%) 95 107 0e 


52% O30. 5s 


= 122% @/ 30.5° 


145.8% @ /-108.4° 


145.8% @ /108.4° 


3 


103 


“Wa 


a 
opt.=  ¢ yaa 70) nee a to ante oan 


| we 
ca’ ane ee. 
: | aa IF 
_ 7 : 7 / is 
» »¢ 1% . wT i ae; . 


See Sed 7 


ae | 


iw a 


=e 
ae) 
i TAS 


= 


2 
ay | ee 


is * @ ie - 
4 y Nts i bi SE a, 


B30191 


